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ABSTRACT 

We present Very Large Array observations of NH 3 (1,1) and (2,2) in a sample of six infrared-dark clouds 
(IRDCs) with distances from 2 to 5 kpc. We find that ammonia serves as an excellent tracer of dense gas in 
IRDCs, showing no evidence of depletion, and the average abundance in these clouds is 8.1 x 10" 7 . Our sample 
consists of four IRDCs with 24 yum embedded protostars and two that appear starless. We calculate the kinetic 
temperature of the gas in IRDCs and find no significant difference between starless and star-forming IRDCs. 
We find that the bulk of the gas is between 8 and 13 K, indicating that any embedded or nearby stars or clusters 
do not affect the gas temperature dramatically. Though IRDCs have temperatures and volume densities on 
par with local star formation regions of lower mass, they consist of much more mass which induces very high 
internal pressures. In order for IRDCs to survive as coherent structures, the internal pressure must be balanced 
by a confining pressure provided by the high concentration of molecular clouds in the spiral arm in which they 
reside. The high molecular concentration and pressure is roughly consistent with gas dynamics of a bar galaxy. 
Subject headings: 



1. INTRODUCTION 

Tracing the evolution of the Galaxy relies heavily on what is 
known of star formation. Here, in the outskirts of the Galac- 
tic disk near the Sun, however, only a part of the picture is 
represented. Low-mass star formation dominates the local 
neighborhood, while the massive stars that dominate the en- 
ergetics in the Galaxy, lie much further away. If we are to 
follow in the footsteps of the successful studies of nearby re- 
gions of star formation, it is essential to isolate and charac- 
terize objects in the evolutionary stages leading the formation 
of a massive star and cluster. This poses a daunting challenge 
because, as we've learned from low mass studies, such ob- 
jects in early stages are cold, dark, and obscured by heavily 
extincted envelopes, and these conditions are only more ex- 
treme in the case of more massive star formation. Thus, we 
must turn from the optical and near infrared studies and look 
longward in wavelength where dust and molecules reveal the 
nature of the early stages of massive star formation. 

Over the past decade, infrared-dark clouds (IRDCs) have 
become the front-runners for the precursors to massive stars 
and clusters (e.g. Carey et al. 2000; Rathborne et al. 2006). 
First, their temperatures, densities, and masses (T < 20 K, 
n > 10 5 cm" 3 , M ~ 10 2 - 10 3 AT©) suit the conditions re- 
quired for massive stars to form (e.g. Ragan et al. 2006; Simon 
et al. 2006). Furthermore, their substructure (Butler & Tan 
2009; Hennemann et al. 2009), including their clump mass 
function (Ragan et al. 2009; Peretto & Fuller 2009), mass- 
radius relationship (Kauffmann et al. 2010) show that IRDCs 
are undergoing fragmentation. Finally, observations of pre- 
and protostellar cores (Henning et al. 2010) within IRDC fila- 
ments, young stellar objects on the IRDC surface (Ragan et al. 
2009), and embedded massive stars in the densest regions (e.g. 
Rathborne et al. 2007; Beuther & Steinacker 2007) show that 
IRDCs are actively forming stars, some of which are massive. 
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This wealth of new insight into the nature of IRDCs has re- 
sulted in large part from the advent of Galactic Plane surveys 
in the mid-infrared (MSX, ISO, Spitzer), where IRDCs appear 
as absorbing structures against the background provided pri- 
marily by polycyclic aeromatic hydrocarbon emission. Fol- 
low up studies using single-dish observations of dust contin- 
uum emission have also been used as a benchmark for the 
gross properties of IRDCs, but such surveys suffer from in- 
sufficient angular resolution to probe the "clump" and "core" 
scales seen with Spitzer. Only a handful of interferometric 
studies have been done (e.g. Hennemann et al. 2009; Beuther 
& Henning 2009), but our capability to conduct such surveys 
will markedly increase with ALMA. In the meantime, there is 
much to be learned about IRDCs from the study of molecular 
line emission. 

The ammonia molecule (NH3) has been widely used to 
characterize the physical conditions in nearby molecular 
clouds (e.g. Wiseman & Ho 1998; Jijina et al. 1999; Li et al. 
2003; Rosolowsky et al. 2008; Friesen et al. 2009). A sym- 
metric top molecule, NH3 is a valuable tool because of its sen- 
sitivity to a broad range of excitation conditions and it probes 
the very densest regions of molecular cloud cores. In partic- 
ular, the lower metastable states of the vibrational inversion 
transitions, (/, K) = (1,1) and (2,2), are excited at the charac- 
teristically low (< 20 K) temperatures of molecular clouds 
and IRDCs. The splitting of the inversion transitions into 
hyperfine structure, due to the interaction between the N nu- 
cleus' electric quadrupole moment and the electric field of the 
electrons, allows for the straightforward calculation of optical 
depth, column density, and gas temperature (see Ho & Townes 
1983, for derivation details). In these nearby clouds, ammonia 
is found to preferentially trace the densest material with tem- 
peratures between 10 and 20 K. It is therefore of interest to 
test whether this correlation holds in the IRDC environment 
where clusters and massive stars are forming. 

Beyond local molecular clouds, Pillai et al. (2006) mapped 
ammonia in nine IRDCs using a single dish telescope with 
resolution ~ 40". They noted that some IRDCs show spatial 
temperature gradients on large scales. However, these obser- 
vations do not resolve the dense cores and clumps that ex- 
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ist on 10 " scales (Ragan et al. 2009). Furthermore, Spitzer 
has revealed the presence of young stellar populations on the 
surface of and within IRDCs, and the possible resulting tem- 
perature variation on these small spatial scales has not been 
examined systematically. 

We present observations towards a sample of six infrared- 
dark clouds from the (Ragan et al. 2009) sample of the NH 3 
(/, K) = (1,1) and (2,2) inversion transitions using combined 
data from the Very Large Array (VLA) and Green Bank Tele- 
scope (GBT). We relate the distribution of ammonia to the 
column density distribution in IRDCs, as well as the absence 
and presence of embedded star formation, from our Spitzer 
study (Ragan et al. 2009). This sample exhibits a mixture of 
apparently starless IRDC clumps and regions of active star 
formation, and we thus relate the star formation properties to 
signatures seen in the lines. Using the unique properties of 
the ammonia transitions as a gas thermometer (Walmsley & 
Ungerechts 1983), we map the kinetic temperature of the gas 
at angular resolution between 4 and 8", corresponding to 0.08 
and 0.16 pc spatial scales at 4 kpc. This study addresses the 
changes in gas temperature on the spatial scales relevant to 
cloud fragmentation. The kinematic characteristics of these 
objects are not analyzed here, but will be addressed in a fu- 
ture paper. 

2. OBSERVATIONS & DATA REDUCTION 

We select a subsample of six IRDCs from the Ragan et al. 
(2009) survey. This sample is biased toward the most compact 
and opaque objects, in hopes that this increases the possibility 
of isolating massive starless (or pre-stellar) cores and clumps. 
The observation details can be found in Table 1. We obtained 
both single dish and array data of the NH 3 (1,1) and (2,2) 
lines. 

2.1. Green Bank Telescope 

We acquired single-dish observations of NH3 (J, K) = (1,1) 
and (2,2) inversion lines using the Robert C. Byrd Green 
Bank Telescope (GBT) from 6 to 15 September 2006. The 
rest frequencies of the (1,1) and (2,2) lines are 23.6944955 
and 23.7226336 GHz (Ho & Townes 1983), respectively. At 
these frequencies, the beam FWHM was approximately 32". 
These upper K-band observations were made in frequency- 
switching mode. The GBT spectrometer back end was config- 
ured to simultaneously observe the two transitions in separate 
spectral windows, using -32,768 spectral channels in each IF 
that were 50 MHz wide. The spectral resolution for the (1,1) 
line was 0.03 km s" 1 , or 1.6 kHz. 

We mapped the regions using a single-pointing grid with 
Nyquist sampling, the parameters of which are summarized 
in Table 1. Integration times were 2-3 minutes per point, de- 
pending on the strength of the line, elevation of the source, 
and thus the time required to obtain sufficient signal-to-noise. 
Pointing corrections were done toward various calibrators de- 
pending on the time of observation (1733-1304, 1814-1853, 
1822-0938, 1831-0949, and 1851+0035) every 45 to 75 min- 
utes, weather-depending, which resulted in corrections of a 
few arcseconds, typically. System temperatures were between 
60 and 100 K, and the elevation of the sources ranged between 
20 and 40 throughout the observations. 

Data were reduced and calibrated using GBTIDL 4 . The 
frequency- switched observations were reduced with the 

4 GBTIDL is an interactive package for the reduction and analysis of spec- 
tral line data taken with the GBT. http://gbtidl.nrao.edu/ 



getf s routine, which retrieves, calibrates, and plots the spec- 
trum. We then applied hanning smoothing and then boxcar 
smoothing over 50 channel bins. Five of the spectral compo- 
nents of the NH3 signature were always spectrally resolved, 
and we performed gaussian fitting on each hyperfine compo- 
nent. The data were then put into a FITS data cube using a 
homemade IDL-based script. 

2.2. Very Large Array 

Observations of the NH3 (1,1) and (2,2) inversion transi- 
tions were undertaken of the sample in the compact D con- 
figuration of the Very Large Array (VLA) over the course 
of three tracks in 2007 April. At this time, the array con- 
sisted of a hybrid of VLA and the updated EVLA receivers, 
and Doppler tracking was unavailable. The observations were 
conducted in fixed-frequency mode. Bandpass and calibration 
were done with observations of 1331+305 and 2253+161. 
Phase calibration was done with observations of 1833-210, 
1832-105, and 1851+005 every -12-15 minutes, and point- 
ing was also performed on these sources every ~ hour. The 
data were flagged and calibrated using NRAO Astronomical 
Image Processing System (AIPS 5 ) according to the method 
described in the AIPS cookbook, taking into account correc- 
tions for the absence of Doppler tracking. 

These K-band observations were made using a four IF 
correlator backend configured with 3.125 MHz bandwidth 
containing 64 channels, yielding a spectral resolution of 
48.8 kHz, or 0.6 km s" 1 . This setup was selected as to ensure 
detection of relatively weak lines in IRDCs, both the (1,1) 
and (2,2) simultaneously, with enough spectral resolution to 
resolve the lines. Unlike the GBT setup, this configuration 
only fits the central three of the five main components of the 
ammonia signature (see Figure 1). We discuss the effects of 
this in Section 3. Table 1 gives the target, pointing and sensi- 
tivity information for each IRDC. 

2.3. Single -dish and Interferometer Data Combination 

The detectable size of structures has an upper limit when 
using an interferometer, which is set by the shortest spacing 
between two antennas in the array. Any flux on scales larger 
than this upper limit will be missing from the final image, a 
phenomenon known as the short-spacing problem. Because 
IRDCs are complex structures, with emission on all scales, it 
is necessary to recover any large scale emission by combining 
the interferometer data with single-dish observations at the 
same frequency. The 100-meter Green Bank Telescope offers 
the perfect complement to our D array observations with the 
Very Large Array. 

Following Friesen et al. (2009), the GBT were first scaled 
and regridded to match the resolution of the VLA data. Using 
the MIRIAD task IMMERGE, the deconvolved interferometer 
images were combined in the Fourier domain with the GBT 
images. IMMERGE uses a tapering function which weights 
short spacings higher than long spacings. Using an annulus 
ranging from 30 to 70 meters, where the high and low resolu- 
tion images overlap, we ensure accurate determination of the 
flux calibration factor. The combined data sets were binned 
in velocity to 0.6 km s" 1 , which reduced the level of random 
errors. This method was applied to the (1,1) and (2,2) datasets 
independently, yielding images of the same spatial and veloc- 
ity resolution as the VLA image. Table 1 gives the final beam 

5 http://www.aips.nrao.edu/ 
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TABLE 1 

Target Summary. IRDC positions, distance, rms and beamsizes of the combined dataset. 



IRDC 




a 


8 


D 


Vfar 


rms 


beam size 






(J2000) 


(J2000) 


(kpc) 


(km s" 1 ) 


(mJy) 


(" x ") 


G005.85- 


-0.23 


17:59:53 


-24:00:10 


3.14 


17.2 


2.8 


7.7 x 6.8 


G009.28- 


-0.15 


18:06:54 


-20:58:51 


4.48 


41.4 


4.8 


8.3 x 6.4 


G009.86- 


-0.04 


18:07:40 


-20:25:25 


2.36 


18.1 


4.3 


8.1 x6.3 


G023.37- 


-0.29 


18:34:51 


-08:38:58 


4.70 


78.5 


2.5 


5.7 x 3.7 


G024.05- 


-0.22 


18:35:52 


-08:00:38 


4.82 


81.4 


4.3 


8.2 x 7.0 


G034.74- 


-0.12 


18:55:14 


+01:33:42 


4.86 


79.1 


6.8 


8.1 x7.0 



G005.85-0.25 peok position 




Fig. 1. — Example spectrum from the central position of G005. 85-0.23 from 
the GBT, VLA, and the combined data set. Overplotted in the dashed green 
line on the combined spectrum is our fit to the central three components of 
the NH 3 (1,1) signature. 

dimensions, as well as other source parameters, for each ob- 
ject. 

3. RESULTS 

The combined data were made in to spectral line maps in 
both the (1,1) and (2,2) transitions in a single VLA field of 
view (~2' at 1.3 cm) for each IRDC. Spectra at each position 
were fit with a customized gaussian minimization routine, and 
the peak flux, centroid velocity, and linewidth were measured. 
As the VLA backend configuration did not fully cover the 
(1,1) spectral signature, when estimating the integrated inten- 
sity, we assume the intensity ratio described in Ho & Townes 
(1983), where the outermost satellite lines account for 0.22 of 
the intensity compared to that of the main line. The total NH3 
(1,1) integrated intensity does not factor into the calculation 
of the column density or temperature discussed in the present 
work. In Table 2 we list the position of the NH3(1,1) emission 
peak, the velocity and linewidth there, and the total ammonia 
abundance and average temperature over the whole mapped 
region. 

3.1. Ammonia as a tracer of high density gas in IRDCs 

Ammonia has long been used as a reliable tracer of high 
density gas in studies of low-mass cores (Jijina et al. 1999; 
Rosolowsky et al. 2008; Friesen et al. 2009) and high-mass 
cores (Harju et al. 1993; Molinari et al. 1996; Wu et al. 2006), 
where other molecules, such as CO, are depleted (Bergin & 
Langer 1997). The initial studies of ammonia in IRDCs (Pillai 
et al. 2006; Wang et al. 2008) show that ammonia is also a 
good tracer in these environments that have the potential to 
form clusters containing massive stars. The ammonia column 
density can be found through analysis of the (1,1) hyperfine 
structure, the derivation of which can be found in Appendix 
A.l. 



In Figure 2, we show ammonia column density contours 
plotted over the 8 jim Spitzer image for each IRDC. The reso- 
lution element is shown in the lower-left corner of each panel. 
In Figure 3, we plot the column density of NH3 as a function 
of H2 column density for each IRDC. The H2 column density 
is calculated from the absorption in the Spitzer IRAC 8 jim 
image 6 (see Ragan et al. 2009). For these analyses, the H2 
column density map was convolved to match the VLA reso- 
lution for each IRDC, thereby lowering the inferred column 
(see Vasyunina et al. 2009). In Figures 3, we do not plot data 
from regions where the NH3 line opacity exceeds 5, as it is 
not used in our abundance fits (see below). 

In most cases we find a one-to-one relationship between 
these two quantities. The slope of the linear relation is an 
estimate of the ammonia abundance, xnh 3 • Our fits only in- 
clude regions in the cloud that have low line optical depths 
(r m (l, 1) ~5). Individual abundance values are given in Ta- 
ble 2, except G009.28-0.15 for which no linear fit was pos- 
sible. In this case, although most points (in black) show a 
one-to-one relation similar to the other IRDCs. We inves- 
tigated where the apparently anomalous emission originated 
and found that "high" values (in green) are confined to a 
southern region and the "low" values are computed for the 
northern tip. Apart from being contiguous, these datapoints 
are not otherwise distinct from the rest of the cloud (i.e. not 
an optical depth effect) or the rest of the sample, therefore we 
do not discard them for the sake of the desired fit. Our restric- 
tion on the line opacity most heavily affects the data fit for 
G023.37-0.29, which excludes all of the central region where 
the NH3 (1,1) lines are saturated. Ragan et al. (2009) also re- 
port the highest values of rg which suggests that N(H2) in the 
central part of G023.37-0.29 is also an underestimate. Thus 
the abundance here is calculated for the outer edges of the 
cloud where both tracers are least hampered by optical depth 
effects. 

The IRDCs G005.85-0.22, G009.86-0.04, G024.05-0.22, 
G034.74-0.12, and the outer regions of G023.37-0.29 show 
good correlation between the column densities of H2 and am- 
monia. The average of these five IRDC abundances is 8.1 x 
10" 7 , plotted in red on each panel. IRDCs with no evidence 
of embedded star formation activity show the tightest relation 
between ammonia and H2 column densities. 

Our abundance value is higher than 4 x 10" 8 found by Pil- 
lai et al. (2006) in a single-dish study of ammonia in IRDCs. 
Though we report similar H2 column densities, we probe 
higher NH3 column densities, namely because our study has 
effective beams at least five times smaller than single-dish 
studies. In a single-dish study, the compact emission is beam- 
diluted which would result in lower NH3 intensity measure- 

6 The H2 column densities fall below the 10 23 cm -2 "saturation limit" 
reported by Peretto & Fuller (2009), indicating that it is a faithful tracer of 
H 2 in these IRDCs. 



4 



Ragan, Bergin, & Wilner 




Fig. 2. — NH3 column density contours plotted over the Spitzer 8//m image of the IRDC. The magenta crosses mark the locations of the emission peaks described 
in Table 2. The other symbols signify the locations of YSOs as follows: Class I (blue square), Class II (green diamond), transition disk objects (cyan circles), 
embedded objects (red diamonds), and unclassified point sources appearing only at 24 /urn (star symbols). Top Left: G005. 85-0.23 N(NH3) contour levels: 2 to 
10 x 10 15 by 1 x 10 15 cm" 2 . Top Middle: G009.28-0.15 N(NH 3 ) contour levels: 3 to 17 x 10 15 by 2 x 10 15 cm" 2 . Top Right: G009.86-0.04 N(NH 3 ) contour 
levels: 2 to 6 x 10 15 by 1 x 10 15 cm" 2 . Bottom Left: G023.37-0.29 N(NH 3 ) contour levels: 3 to 18 x 10 15 by 5 x 10 15 cm" 2 . Bottom Middle: G024.05-0.22 
N(NH 3 ) contour levels: 2 to 10 x 10 15 by 1 x 10 15 cm" 2 . Bottom Right: G034.74-0.12 N(NH 3 ) contour levels: 8 to 30 x 10 15 by 2 x 10 15 cm" 2 . 



TABLE 2 

Peak properties, abundances, and mean temperature of IRDCs 
Peak NH 3 f Tdv position 



IRDC 




RA 


DEC 


V 


dv 


r m (l,l) 


XNH 3 


<T k > 






(J2000) 


(J2000) 


(km s" 1 ) 


(km s" 1 ) 




(10- 7 ) 


(K) 


G005.85- 


-0.23 


17:59:51.4 


-24:01:10 


17.4 


1.4 


4.9 


4.4 


11.0+1.0 


G009.28- 


-0.15 


18:06:50.8 


-21:00:25 


41.0 


1.5 


4.0 


a 


11.2+1.3 


G009.86- 


-0.04 


18:07:35.1 


-20:26:09 


17.8 


1.4 


3.3 


5.2 


10.3+1.1 


G023.37- 


-0.29 


18:34:54.1 


-08:38:21 


78.5 


3.7 b 


b 


5.4 


8.8±3.3 


G024.05- 


-0.22 


18:35:54.4 


-07:59:51 


81.7 


2.0 


2.6 


6.7 


10.1 + 1.4 


G034.74- 


-0.12 


18:55:09.5 


+01:33:14 


78.0 


2.3 


6.1 


19. 


8.5±0.8 



a Linear fit did not converge. 

b Line is saturated and optically thick at emission peak. Abundance fit in 
cloud edges. 

ments than what we find here with the VLA. 

3.2. Temperature maps 

The temperature in cold molecular gas is attainable with 
knowledge of the optical depth of the line and the ratio of 
NH3 (1,1) and (2,2) central line intensities. The derivation of 
the kinetic temperature is derived in Appendix A. 2. 

This sample contains objects both with and without signa- 
tures of star formation. In Figure 4, we present the kinetic 
temperature maps in all six IRDCs with contours showing the 
intensity of NH3 (1,1) and symbols indicating the locations of 



young stars. The class II sources (the most evolved YSOs dis- 
played) are the most widely distributed population and do not 
follow the dense gas traced by ammonia. As such, they are 
unlikely to have a strong impact of the gas compared to the 
embedded protostars. The temperature is only calculated for 
positions with detections three times the rms noise. Since the 
(2,2) emission is typically much weaker than the (1,1) emis- 
sion, the temperature map boundaries are limited by the (2,2) 
signal-to-noise. Kinetic temperatures range from 8 to 13 K. In 
Figure 5, we present the kinetic temperature as a function of 





Fig. 3. — NH3 column density versus the H2 column density (10% error bars shown). The red line shows the trend of the linear fit of the abundance, xnh 3 - 
N(NH3)/N(H2). The average abundance in this sample is 8.1 x 10 -7 . The fits were only performed in unsaturated regions of moderate total optical depth in the 
(1,1) line (t(1, 1, m) < 5). For G009.28-0.15, a linear fit did not converge. We show in green and cyan the values corresponding to the regions circled in green 
and cyan in Figure 2. These regions appear to have anomalously high or low NH3 emission (respectively), and do not keep to the one-to-one relation the other 
IRDCs exhibit. G023. 37-0.29 is only fit for the un-saturated region around the cloud edges. 



the molecular hydrogen column density (N(H2)) derived from 
8 jim absorption 7 (Ragan et al. 2009). Table 2 lists the average 
temperature associated with each IRDC. 

IRDCs G005.85-0.23 and G024.05-0.22 lack any 24 jim 
point sources coincident with the gas traced by NH 3 . We infer 
that these IRDCs contain no embedded star formation, though 
there are stars projected in the vicinity (shown on Figure 4). 
In both cases, the NH3 intensity peaks at the same position 
of both the 8 and 24 jim absorption peaks. The average tem- 
perature profiles (Figure 5) in these two starless IRDCs are 
the flattest of the entire sample (<T^-„ > varies by ~1 K in 
both cases), which is consistent with these objects being least 
influenced by radiation sources. 

The four remaining IRDCs in the sample (G009.28-0.15, 
G009.86-0.04, G023.37-0.29, and G034.74-0.12) coincide 
with at least one 24 jim point source, indicating that these 
objects likely host embedded protostars. The "northern" re- 
gion of G009.28-0.15 (shown in cyan in Figures 2 and 3), 

7 Were IRDCs round, NXH2) could be used as proxy for radius, thus pro- 
viding a radial temperature profile. However, these objects are not spherical 
and tend to be filamentary in nature, so regions of high optical depth are those 
most akin to the "center" of a dark core or clump, and low optical depth re- 
gions are those nearer to edges or near radiation sources such as embedded 
stars. 



which coincides with two 24 jim point sources, appears to 
have a slightly elevated temperature compared to the "south- 
ern" region (shown in green in Figures 2 and 3), which ap- 
pears infrared-dark. While these regions are each spatially 
coherent, the absolute difference in the mean temperature is 
within the temperature measurement error shown in Figure 5 8 . 

There are six 24 /mi point sources in G009. 86-0.04, five 
of which are concentrated in the eastern portion of the cloud 
where it is not possible to calculate the temperature because of 
very weak (2,2) emission. The kinetic temperature at the inte- 
grated intensity peak is 9 K, and the point source in the to the 
southwest has an elevated temperature of 13 K. We see the 
strongest temperature gradient in this source (see Figure 5), 
with the lowest temperatures tending toward the highest opti- 
cal depths. 

IRDC G023. 37-0.29 exhibits the strongest NH 3 emission in 
the entire sample, and we also estimate the highest abundance 
in this object. The lines near the center of this object are satu- 
rated, indicating that the emission is very optically thick. The 
temperature then is likely not representative of the conditions 
near to the 24 jim source near the center of the object, which 

8 In order to minimize the random errors and retain the largest possible 
spatial expanse of reliable temperature measurements, we binned the spectra 
to 0.62 km s _1 resolution. 
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could explain why we do not see any elevation in kinetic tem- 
perature like we do in the other IRDCs in the sample. It is 
likely that this object would potentially exhibit strong emis- 
sion in higher /, K lines and, therefore, be at a higher temper- 
ature not probable by the present observations. 

In the IRDC G034.74-0.12, emission is widespread in both 
the (1,1) and (2,2) transitions. We find the highest ammonia 
column densities in this object (~ 3x 10 16 cm -2 . This map 
is also the noisiest, and the (2,2) emission map only passed 
our 3-sigma detection limit test at the two strong peaks also in 
(1,1) emission. If the noise requirement for the (2,2) emission 
is relaxed, the resulting temperature across the source is very 
uniform, between 8 and 9 K. Four 24 yum point sources are de- 
tected on the IRDC, and the IRAS point source 18526+0130 
(Lboi ~ 10 4 L o Robitaille et al. (2007)) is 30" north of the 
NH 3 map boundary, corresponding to 0.7 pc at the distance to 
this IRDC, though no direct influence of the IRAS source is 
evident in our temperature measurements. 

4. DISCUSSION 

4.1. Column density structure 

Throughout studies of star formation regions with molec- 
ular emission, ammonia has been one of the best tracers of 
cold dense gas, showing no compelling evidence of depletion 
up to very high densities, n > 10 5 cm -3 (Bergin & Langer 
1997) . The same holds true in IRDCs, though in regions of 
very high H2 column density, the ammonia line can become 
optically thick. We find no evidence of depletion of ammonia 
in IRDCs, thus it can serve as a reliable tracer of kinematics 
throughout the bulk of the gas which is not yet forming stars. 
We will, in the following sections, discuss the limitations of 
the (1,1) and (2,2) lines in probing warm gas in compact re- 
gions near embedded protostars. 

We note that our estimated NH3 abundances are high when 
compared to other cold objects. For example, Tafalla et al. 
(2004) find central abundances on the order of ~ 2 x 10" 8 in a 
few cold low mass cores in Taurus, while Harju et al. (1993) 
find values from 1 - 5 x 10" 8 in a sample of isolated cores in 
Orion and Cepheus. Interestingly, Tafalla et al. (2004) find ev- 
idence for an increase in the NH3 abundance in the denser gas. 
They postulate that perhaps greater CO depletion might lead 
to enhancements of the ammonia abundance via gas phase re- 
actions. Our sources likely have greater densities and perhaps 
that induces greater formation of ammonia. However, this 
question needs to be addressed with more complete modeling 
spanning a wide range of densities. 

4.2. IRDC temperatures 

We find that the temperatures in this sample of IRDCs, de- 
spite spanning from quiescent to active phases, are largely 
uniform in the bulk of the cloud, ranging only from 8 to 13 K. 
These results are in general agreement with the single-dish 
measurements of temperature in IRDCs by Pillai et al. (2006), 
and also the temperatures found in Perseus (Rosolowsky et al. 
2008) and Ophiuchus (Friesen et al. 2009). 

Ammonia studies of temperature structure in massive star 
forming complexes have mainly focused on Orion. Wise- 
man & Ho (1998) map the filamentary region surrounding the 
Orion- KL massive star forming core complex. One source 
within this complex is IRc2, which drives a large bipolar 
outflow. The surrounding molecular ridge, OMC-1, is com- 
prised of several finger-like filaments in which they find the 
kinetic temperature ranges between 15 and 30 K. In the cen- 
tral core around Orion-KL, they find temperatures near 80 K, 



and near the edges of the filaments, the temperature is ele- 
vated to -40 K. The heating is detectable up to 0.3 pc away 
from the central core, and is attributed to radiation and winds 
from Orion-KL sources and shock interactions near the core 
where the highest temperatures are found. If such luminous 
sources were present in IRDCs, we would see greater temper- 
ature variation. The low, unvarying temperature of the bulk of 
the IRDC gas indicate that IRDCs are in an earlier evolution- 
ary stage than Orion. 

Li et al. (2003) select quiescent cores far away from the 
Orion-KL region, a site of active massive star formation. The 
authors observe a similar temperature trend as the one pre- 
sented here: the lowest temperatures are found in the cloud 
interiors and warmer gas is near the edges. They attribute this 
to external heating from the elevated radiation field due to 
nearby active star formation region in the Trapezium, despite 
the large amount of extinction surrounding the region probed 
by ammonia. These quiescent Orion cores can be most favor- 
ably compared to the IRDCs G005.85-0.23 and G024.05-0.22 
in our sample. They show no 24 yum point sources coincident 
with the dense gas. We see very little heating in these ob- 
jects, and by comparing the colored region (where there is a 
NH 3 (1,1) and (2,2) detection) and the contours (NH 3 (1,1) 
intensity), the region where temperature calculation is possi- 
ble is already heavily embedded in dense gas, and therefore 
shielded from the interstellar radiation field. 

One high-resolution follow-up study of IRDC 
G028.34+0.06 Wang et al. (2008) showed a significant 
rise in rotation temperature at the sites of 24 yum point 
sources. In this case, the more evolved, warmer region 
contains a well characterized high-mass protostellar object 
and IRAS source (L ~ 10 3 L o ) and has a peak rotation 
temperature of 30 K, much higher than detected in any IRDC 
in this study. In the single-dish study of this object, however, 
the rotation temperatures across the IRDC were found to be 
between 13 and 16 K. Thus, Wang et al. (2008) demonstrate 
that without the high-resolution of an interferometer, these 
compact features are diluted by the emission from the cold, 
quiescent bulk of the IRDC. 

Our sample of six IRDCs was selected to be at the earliest 
evolutionary stage, in regions devoid of maser emission, radio 
emission, and IRAS sources (Ragan et al. 2006). Thus, the 
lower average kinetic temperature (~11 K) is consistent with 
this selection, and this sample of IRDCs are at an earlier stage 
in their evolution compared to G028.34+0.06. Furthermore, 
we have no instances of such elevated temperatures such as 
that seen in Wang et al. (2008), indicating that nowhere in 
these IRDCs is there massive star formation at such an ad- 
vanced stage, if it exists at all. 

The absence of temperature structure does not preclude the 
presence of star formation. Indeed there exists a surface popu- 
lation of young (low-mass) stars in these IRDCs (Ragan et al. 
2009), however, these appear to have a minimal impact on the 
gas traced with these observations. In the cases where we de- 
tect embedded star formation, the effect on the temperature is 
localized and does not influence the widespread dense gas in 
IRDCs 9 

4.3. Fragmentation scale constraints 

9 IRDCs show high extinction levels, with typical Ak values in the 1 - 
3 range (or A24 from 0.5 to 1.5 Flaherty et al. 2007) for associated YSOs 
(Ragan et al. 2009) and higher still for more deeply embedded sources not 
detected by Spitzer. For example, assuming A24=l, or A^=2, at 5 kpc, up to 
a 100 L Q star could fall below our detection limit. 




Fig. 4. — Kinetic Temperature map of IRDCs. Overplotted are contours of the NH3 (1,1) integrated intensity, beginning at 0.2 Jy beam 1 km s 1 and increasing 
in steps of 0.1 Jy beam -1 km s -1 . The 0.2 pc linear scale is indicated on each panel. Star symbols indicate the location of 24 jum point sources. Other symbols are 
as follows: red diamonds are embedded objects, green diamonds are Class II objects, and blue squares are Class I objects, all as classified in Ragan et al. (2009). 
Only one known IRAS source is present in the fields of the sample: IRAS 18526+0130 which is labeled in the lower right panel near G034.74-0.12. 




The temperatures determined in this study can place con- 
straints on how fragmentation proceeds in IRDCs. The tra- 
ditional characteristic fragmentation scale (which assumes 
isothermality, spherical symmetry, and uniform density) is 
given by the Jeans length: 



(1) 



where c^ is the isothermal sound speed (0.18 km s" 1 for 10 K 
gas). For an average density of 10 5 cm -3 (p ~ 4 x 10" 19 
g cm -3 ), the Jeans length is 0.07 pc, which is below the 
effective physical resolution of our VLA observations (0.1 
to 0.2 pc). The corresponding Jeans mass (M/ = pA 3 j) is 
about 2 M Q . While the detectable star formation activity 
within these IRDCs is limited to a few objects, in the cases 
where there are multiple objects of class I or earlier (G009. 28- 
0.15, G009.86-0.04, G034.74-0.12), the average separation 
between protostars is 0.4 pc, within a factor of 5 of the Jeans 
length scale. This results is similar to what Teixeira et al. 
(2006) found in a filamentary star formation region (with a 
factor of three of Aj), which is roughly consistent with ther- 
mal fragmentation. Given that IRDCs are a several kiloparsec 
distances, hampering our sensitivity to faint stars, our scale is 
an upper limit. 

In reality, the material in IRDCs is much more clumped and 
filamentary than a uniform sphere. An alternative quantifica- 
tion of such a cloud is the critical fragmentation mass-per- 



unit-length scale in a filament (e.g. Inutsuka & Miyama 1992; 
Andre et al. 2010), which gives 15 M© pc -1 for 10 K gas (not 
taking into consideration any magnetic field support). The 
mass-per-unit-length in these IRDC filaments are typically a 
few hundred solar masses per parsec, well in excess of the 
critical value. This would indicate that these objects are un- 
dergoing collapse and actively forming stars. The dynamical 
state of this IRDC sample will be analyzed in detail in a forth- 
coming publication (Ragan et al., in prep.). 

4.4. The high pressure environment of IRDCs 

In Ragan et al. (2009) we discussed the overall mass distri- 
bution of IRDCs in comparison to the Orion Molecular Cloud, 
which is the local example of massive and clustered star for- 
mation. From our absorption study we can estimate a gas 
mass surface density for IRDCs, Zirdc, between\000and5000 
M Q pc -2 , while in Orion the mass surface density is ^omc ~ 
30 M pc -2 . Thus IRDCs are not distant Orions, but rather 
something fundamentally different. Here we examine what 
role the environment of IRDCs might play in supporting such 
a large concentration of material where the prevalence of star 
formation appears relatively low. 

Jackson et al. (2008) has examined the spatial distribution 
of IRDCs which were previously identified with the Molec- 
ular Ring that contains over 70% of the galactic molecular 
mass. They argue that the distribution is best fit by a model 
of the Milky Way with two major spiral arms emerging from 
a central bar. The first quadrant IRDCs lie preferentially in 
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Fig. 5. — Kinetic temperature (Tkin) as a function of H2 column density (N(H2)) calculated from 8 //m absorption. The green points indicate the weighted 
average trend of the data. 



the Scutum-Centarus arm. This is in agreement with other re- 
sults probing the dust structure in the galaxy (e.g. Drimmel & 
Spergel 2001 ; Benjamin et al. 2005). 

Theoretical models of gas dynamics in barred galaxies sug- 
gest that gas flow produces arms of gas/dust density enhance- 
ment that extend from the leading side of the bar (Athanas- 
soula 1992). These enhancements are the result shocks and 
other dynamical processes which results in the formation of 
molecular clouds from the atomic medium (Combes & Gerin 
1985; Bergin et al. 2004). The IRDC surface density corre- 
sponds to an average density of < nu 2 >~ 10 5 cm -3 . Given a 
typical uniformly low temperature of ~ 10 K this corresponds 
to a internal thermal pressure of V int jk ~ few xlO 6 cm -3 K 
well in excess of that in Orion, where the gas temperature 
is typically 20 K, or the general ISM where ¥ ism Ik = 2240 
cm -3 K (Jenkins & Tripp 2001). These high pressures require 
a strong dynamical interaction to gather ISM gas (c.f. Bergin 
et al. 2004), enabling gravitational collapse to commence, and 
fragmentation to proceed. This interaction may be related to 
gas flow induced by the central Galactic bar. 

5. CONCLUSIONS 

Ammonia is an excellent tracer of the dense gas in IRDCs. 
To the resolution of our VLA observations, we find a roughly 
linear relationship between ammonia and molecular hydrogen 
column density up to N(H2) ~ 3 x 10 22 cm -2 as traced by 
mid-infrared extinction. We approximate the abundance of 



ammonia in IRDCs to be 8.1 x 10 7 . We see no evidence of 
ammonia depletion in the bulk of the gas in IRDCs. 

Our NH 3 (1,1) and (2,2) observations show that IRDCs 
have uniform temperatures between 8 to 13 K, varying only 
by 2 to 4 K within a given cloud. We see no strong evidence 
for external heating of the gas in IRDCs as probed by am- 
monia. IRDCs are embedded in large complexes of gas that 
shield it from the radiation of local stars. Embedded proto- 
stars, appearing as 24 jim point sources, and nearby young 
stars do not affect the temperature of the bulk gas in IRDCs. 
These IRDCs are drawn from a sample of objects with no 
known signposts of massive star formation (e.g. masers, 
IRAS sources), and the temperature structure also supports 
the absence of such activity. At the same time, the limitation 
of the NH3 (1,1) and (2,2) thermometer in terms of excitation 
compromises the sensitivity to the warm gas that may reside 
on the smallest scales. Observations of more highly excited 
lines would be a useful test of the existence of a warmer, more 
compact core near the embedded sources. 

The gas in IRDCs is undergoing collapse and fragmenta- 
tion, and they will ultimately form clusters containing massive 
stars. Their physical properties (e.g. temperature, density) 
are similar to gas in local star forming regions, but their high- 
pressure environment allows for there to be much more mass 
concentrated in a small area. IRDCs primarily reside in the 
Scutum-Centarus arm as has been suggested in the literature 
and the high concentration of dense high-pressure molecular 
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APPENDIX 
DERIVATION OF PHYSICAL QUANTITIES 
Column Density 

We calculate the NH 3 (1,1) column density following Rohlfs & Wilson (2004), taking into account both parity states of the 
(1,1) energy level (Friesen et al. 2009). 



2 1 8i . p -hv /(kT ex ) 



8ttv^ 1 



where gi and g2 are the statistical weights for the upper and lower states of the NH 3 (1,1) transition (equal in the case of our 
two-level approximation), A(l, 1) is the Einstein spontaneous emission coefficient, 1.68 xl0~ 7 s _1 , and J r(v)dv = ^7ji<j v v§t\Ic 
with vo the rest frequency, T\ the total line opacity, and cr v the dispersion (Rosolowsky et al. 2008). T ex is the gas excitation 
temperature, which can be found by fitting the observed spectrum, 7^ y , for a given frequency 

T a v (J, K) = n MB ®(J(T e M K)) - J(T bg ))(l - exp(-r v (l K))) (A2) 

where t]mb is the main beam efficiency (0.81), O is the filling factor (assumed here to be 1), J(T) is the Rayleigh- Jeans temper- 
ature, and Tb g is the background temperature, 2.73 K. Finally, the total column density of ammonia is determined via use of the 
partition function, Z: 



Z = ^(27+ l)S(J) exp 



-h[BJ(J+ 1) + (C-B)J 2 ] 



kT k 



(A3) 



where B and C are the rotational constants of NH3, 298192 and 186695 MHz respectively 10 , and S(J) accounts for the nuclear 
spin statistics of the system. The total ammonia column density is then 

N(NH 3 ) = Af(l,l)x — (A4) 

in which the Z in the numerator is summed over J=0, 1 and 2, and Z\ in the denominator is evaluated only at J=l. At the 
low temperatures of these IRDCs, we assume these metastable states represent a good approximation to the column density of 
ammonia. 

Temperature 

The rotation temperature derived from NH3 inversion lines depends on the optical depth and the ratio of the main lines in the 
(/, K) = (1,1) and (2,2) signatures. In turn, the optical depth of the NH3 (1,1) main line (r m (1,1)) can be measured by comparing 
the ratio of the main line intensity to the satellite line intensity, as described in Ho & Townes (1983). 

AT a {J,K,m) _ 1 - e -raK,m) 
AT a (J, K,s) ~ 1- e~ ar ^ K ^ 

where the left hand side is the ratio of the intensities of the main and satellite lines in the NH 3 (1,1) signature, and a is the ratio 
of the intensities compared with the main line intensity. With measurements of AT a for the main line and the two inner satellite 
lines (for which a = 0.28), one can solve for r numerically. The total optical depth of the (1,1) line, t\ is twice the main line 
optical depth (Li et al. 2003). The rotational temperature (T Rot ) characterizes the distribution of molecules in these two energy 
levels, which are separated by 41.5 K in energy, denoted as 7q. 



TRot = -To In 



-0.282 
-In 



T W (1,1) 



Ar a (2,2,^) 
Ar fl (l,l,m) 



(A6) 



Since there are no radiative transitions between the (2,2) and (1,1) levels, their population ratio, represented by T Rot , probes the 
importance of collisions, and thus can be used to estimate T&. This assumption is robust in environments such as pre-stellar cores 
and IRDCs, where temperatures are typically 20 K and lower (Walmsley & Ungerechts 1983). At such temperatures, only the 
(1,1), (2,2), and (2,1) states are important, thus Swift et al. (2005) showed that the relationship between the kinetic and rotational 
temperatures is as follows: 

10 http://spec.jpl.nasa.gov/ 
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f Rot 



\ + l±ln(\+0.6e- l51ITk ) 



(A7) 



Tafalla et al. (2004) perform Monte Carlo models including the six lowest levels of para-NH 3 , varying the temperature between 
5 and 20 K and self-consistently calculating the (1,1), (2,2), and (2,1) level populations. They find that the approximation 



[ Rot 



1 - ^ln[l + l.le- 15 -V T **] 



(A8) 



is accurate to 5% or better in the 5 to 20 K temperature range. Recent work by Maret et al. (2009) shows that this calibration of 
the ammonia thermometer holds with this precision taking into account the most up-to-date collision rates in the literature. Thus, 
the two-level approximation for the ammonia thermometer described above is robust. 
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